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Abstract 

00 

£Nj . Discrete models usually represent approximations to continuum physics. Cylindrical con- 

sistency provides a framework in which discretizations mirror exactly the continuum limit. 
1 , Being a standard tool for the kinematics of loop quantum gravity we propose a coarse grain- 

O -1 ing procedure that aims at constructing a cylindrically consistent dynamics in the form of 

^' ( , transition amplitudes and Hamilton's principal functions. The coarse graining procedure, 

bJTJ' which is motivated by tensor network renormalization methods, provides a systematic ap- 

proximation scheme towards this end. 

A crucial role in this coarse graining scheme is played by embedding maps that allow 
the interpretation of discrete boundary data as continuum configurations. These embedding 
C ■ ' maps should be selected according to the dynamics of the system, as a choice of embedding 

maps will determine a truncation of the renormalization flow. 

^ 1 Introduction 
O 

(N 

Discretizations have become a leading tool for the construction of quantum gravity models [I], 
for instance in loop quantum gravity [2j [3] , spin foams [H [5] or Regge gravity [6] . Indeed discrete 
models provide a very effective method to access non-perturbative physics. 

Given a discretization, the question arises how to relate observables from the discrete model 
to observables in the corresponding continuum limit (should it exist). Usually this question 
is handled by considering a family of discretizations, in which we can consider the refinement 
limit. Observables computed from a given discretizations will in general be approximations to 
corresponding observables in the continuum limit. 

An alternative view is provided by the concept of cylindrical consistency, which is a key 
technique in the construction of the continuum Hilbert space of loop quantum gravity [7]. The 
idea here is that discretizations are not seen as approximation but as a selection of a certain 
set of degrees of freedom, for which the discrete model should give the same predictions as the 
continuum model. More precisely discretizations come with a (partial) ordering into finer and 
coarser and embedding mapg^ of coarser into finer discretizations. Cylindrical consistency de- 
mands that observables which can be represented in a given discretization should not depend on 
the choice of finer discretization into which this given discretization is embedded. In other words, 
a prediction drawn from a given discretization should be also valid for any refined discretization 
or the continuum limit. 



1 More precisely the configuration spaces or Hilbert spaces associated to coarser discretizations are embedded 
into those for finer discretizations. 
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So far this concept has been successfully applied to the construction of the kinematical 
Hilbert space [2], including observables and the Hamiltonian constraints [8]. The question 
arises whether we can implement this concept for the dynamics, for instance in the form of a path 
integral (measure) as suggested in [9] see also [10] for related ideas. This would be equivalent 
to a cylindrical consistent construction of transition amplitudes (or the physical Hilbert space). 
On the classical level we would ask for cylindrically consistent family of Hamilton's principal 
functions which are the classical analogues of transition amplitudes [11] . 

In this way cylindrical consistency not only provides the continuum limit of a given theory 
it also allows via the embedding maps a precise interpretation of discrete configurations (states) 
as elements of the continuum configuration (Hilbert) space. 

The concept of a cylindrically consistent dynamics is actually very similar to the concept of 
perfect actions [12] or perfect discretizations. These were suggested for gravitational systems in 
|13[ [14"] and constructed for a number of examples in |14[ \TE[ [TBI [T7] to address the problem of 
broken diffeomorphism symmetry in generic discretizations of gravity [18} \19\ . 

Indeed recent work [14} \W\ \2U[ [21] has shown that demanding diffeomorphism symmetry 
is a very strong requirement: there are many examples where it is equivalent to demanding 
triangulation (or more generally discretization) independence!! This includes coarse graining 
the triangulation or making one region much denser triangulated than another region. Such a 
notion implies that the discrete theory should correctly reproduce physics on all scales, not only 
on scales large compared to the discretization scale [21] . In other words such a discretization has 
to mirror continuum theory exactly, which is the defining property of a perfect discretization. 

A way to construct such perfect discretizations is as fixed points of a (real space) Wilsonian 
renormalization flow. Such a flow generates non-local couplings (in non-topological theories), see 
for instance |22[ 115], These non-local couplings are crucial in order to obtain a non-topological 
theory which is triangulation independent. 

The non-local couplings make the theory also computationally very challenging. As one 
basically attempts to solve the dynamics of the theory, a systematic framework that would im- 
prove the discretization in steps, and also allows a control on the amount of non-local couplings, 
is advisable. In this work we will propose such a framework. As a key input we introduce 
boundaries, the fineness of the boundary data will determine the quality of the approximation. 
This already indicates that we will have a partial ordering into finer and coarser boundaries and 
corresponding embedding maps. This will allow us to formulate the aim of achieving a perfect 
discretization as finding a cylindrically consistent effective action or Hamilton's principal func- 
tion/ path integral measure. As we will see the choice of embedding maps will be essential and 
should be determined by the dynamics of the system. 

Another view point of how this proposal differs from the usual understanding of perfect 
discretizations, is that in the latter we keep the building blocks of the discretization simple but 
allow for non-local couplings between these building blocks. Here we will introduce more and 
more complicated building blocks with finer and finer boundaries. The couplings between these 
building blocks will still be local, i.e. only variables associated to neighbouring building blocks 
will be coupled. Note however that the number of such variables will be growing and that the 
couplings in-between the variables associated to one building block will be non-local (if we see 
this building block as an extended object). In fact this shift of non-localities from in-between 
building blocks to introducing more and more variables to one building block is inspired by 
tensor network renormalization [231 l24"[ [25] . There the path integral amplitude (together with 

2 This notion of diffeomorphism invariance implies that the positions of vertices in a given solution are gauge 
parameters and can be chosen freely. But if we are allowed to move vertices we can produce very densely 
triangulated regions and very coarse triangulated regions. In the extreme case one can move vertices on top of 
each other thus even reducing the number of vertices in a given discretization, or the number of simplices in a 
given triangulation. 
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the measure) associated to a certain region or building block is encoded into a tensor. Summing 
over (bulk) variables correspond to the contraction of tensors into a new tensor, which typically 
is of higher rank, representing more complicated boundary data. The key question is then how 
to approximate a tensor of higher rank with a tensor of the original rank, that is one has to 
choose a truncation of the renormalization flow. Here we will put forward the question whether 
and how we could use the embedding maps, which underly cylindrical consistency, for such a 
truncation. This relation of embedding maps to truncations of the renormalization flow makes 
it obvious that the embedding maps should be determined according to the dynamics of the 
system. 

In the next sections [2] and [3] we will introduce the concept with the help of an example, the 
scalar field with and without potential in two dimensions. In the following sections and [5] we 
will consider variations of the coarse graining scheme and compare two different choices for the 
embedding maps. The coarse graining scheme will be formalized in section [6j which will also 
shortly review the concept of cylindrical consistency. In section [JJ we will extend the procedure 
to the quantum case and discover that the choice of embedding maps is even more important in 
this case. We end in section [8] with a summary an outlook. 

2 The massless scalar field 

As an instructive example consider the massless scalar field on a triangulation of flat Euclidean 
two-dimensional space. Ond discretization for the scalar field action on a triangulation is given 
by 

S = ^2 Sa = ^2^2 C0t ( a e)(</>s(e) - 4>t(e)) 2 , (2-1) 
A A ecA 

where a g is the angle at the vertex opposite the edge e, and s(e),t(e) denote the source and 
target vertex of e respectively. Note that in two dimensions the action for the massless scalar 
field is scale invariant. 

The action (|2. 1[) is given as a sum over actions associated to the triangles. This additivity of 
the action will be essential for what follows. As we will see the action associated to one triangle 
can be interpreted as Hamilton's principal function for this triangle. 

The dynamics defined by the action (|2.ip has the remarkable property to be invariant under 
the so-called 1-3 Pachner move. This move divides one triangle into three by introducing an 
inner vertex, see figure [TJ (This has to be done such that the deficit angle at the inner vertex is 
zero, i.e. after the subdivision the outer triangle is still flat.) Hamilton's principal function H, 
the action for the complex of three triangles, evaluated on a solution with given boundary data, 
will then agree with the action for one triangle (which we identify with Hamilton's principal 
function for this triangle) 

Haifafafa) = Hxfatfafa) . (2.2) 

Moreover the left hand side does not depend on where we place the inner vertex, that is the 
result is invariant under vertex displacement, which can be identified with a discrete notion of 
diffeomorphism symmetry [271 (HI IS] • 

The reason for this invariance is the simple solution for the field on the inner vertex given 

by 

(/>o(xi,x 2 ) = Xt(j)l + X 2 (f>2 + X S (xi,X2)(f>3 ( 2 -3) 
3 There are a number of discretizations for the 2D Laplacian, but there is not a naturally preferred one [26] . 
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Figure 1: Subdivisions of a triangle. 

where the Xj are barycentric coordinates of the triangle, i.e give the position of the inner vertex as 
a weighted average of the outer vertex positions, with X1+X2+X3 = 1, so xa(xi, X2) = 1 — x\ — x-i- 
This solution is a linear function in the coordinates and hence is also a solution of the continuum 
Laplace operator. For such a linear solution we can exactly give three boundary data which define 
the derivatives in the two coordinate directions. 

We can define this solution (|2.3p not only for the coordinate values of the actual inner vertex 
- we can also take this as a continuum solution. This is the reason why Hamilton's principal 
function will be also invariant under further subdivisions of the triangle, including subdivisions 
of the boundary, as long as the boundary data are chosen to be 'edge wise' linear, i.e. 

cpk = xi<f>i + x m <j) m . (2.4) 

(Note that the set of coordinates with Xi = describes the line through the edge opposite the 
vertex i.) The boundary data we impose are the ones following from the dynamics, as these 
data are induced on the boundary by the solution f)2.3f) . 

We conclude from this discussion that edge wise or piece wise linear boundary fields seem to 
be dynamically preferred. 

A two-dimensional triangulation can be changed into any other topologically equivalent 
triangulations by just applying 3 — 1, 1 — 3 and 2 — 2 Pachner moves. Hence to obtain full 
triangulation independence for our discrete theory we would just need invariance under the 2-2 
move, which changes the way how two triangles are glued into a quadrangle. However the action 
(|2,1|) is in general not invariant under the 2 — 2 move (as the actions before and after the move 
display couplings along the two different diagonals). Additionally demanding invariance under 
the 2 — 2 Pachner move and insisting on simple building blocks will lead to a non- locaS action 

[SUES]. 

We will therefore employ a different strategy, allowing for more and more complicated build- 
ing blocks. The discussion is easier on a regular lattice, so we start with actions on squares 
(which for simplicity we take to be also geometrical squares, i.e. equilateral). The justification 
for this choice is that our aim is to give a construction principle for discrete theories, such that 
the dynamics described does not depend on which choice of discretization or lattice one starts 
with, so we can choose a simple, regular lattice. 

Hence we consider a lattice built out of squares. The action for one square can be obtained 
from the action for a triangular subdivision (12. ip . For a regular lattice we have (the lattice 
constant drops out as the action is scale invariant) 

s = E 5 ° = E E ■ ( 2 - 5 ) 

□ □ i=l,..A 

where i = 1, . . . , 4 is a cyclic ordering of the vertices associated to the square □, so addition of 
4 The couplings will decay exponentially with the distance |22| . 
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Figure 2: The boundary fields on a square. To find the fixed point action with eight boundary fields 
{^itJij} one would have to set i.e. 7112 = as it encodes the deviation from a piecewise linear field. 

indices is to be understood mod 4. The coupling along the diagonal drops out as cot (a) = for 
a orthogonal. 

Now lets subdivide the square into four smaller squares. We obtain eight boundary fields 
and the field on the inner vertex cpo, with notation explained in figure [2j The eight boundary 
fields are given by the four fields on the corners of the composite square, which we denote by 
4>i,i = 1, ... 4, and the boundary fields on the inner vertices of the subdivided edges which we 
denoted by </>y for the vertices between the corners i and j. 

We solve the equation of motion for the inner vertex and insert this solution into the action 
for the four squares, which yields Hamilton's function depending on the eight boundary fields. 
Now we have to compare this to the original action, which is just a function of the four boundary 
fields. Here we use the crucial insight from the previous discussion, which showed that piecewise 
linear fields define a special class of boundary data. So we define a new action as the Hamilton 
function, which we just computed, evaluated on piecewise linear boundary fields, that is we set 
fiij = ^(0* + <A?) f° r the fields on the midpoints of the edges. The resulting action 

S' n = (01-02) 2 + (02-03) 2 + (^3-04) 2 + (04-0l) 2 -a(01-02 + 03-04) 2 

= S u - «(0i - <t>2 + 4>3 - <M 2 , (2-6) 

where a = s for this first step, just differs by a total square from the initial action, introducing 
in particular couplings along the diagonal. 

This process can be repeated inducing a flow which only affects the value of a. The fixed point 
a* = I can be readily found. This fixed point action enjoys an invariance under subdivisions 
of a certain type under the condition that the additional boundary data that result are fixed as 
(specific linear) functions of the initial boundary data. Later we will discuss in which sense this 
action already encodes the continuum dynamics for very special boundary data. 

Moreover this fixed point is unique in the following sense. Parameterizing all quadratic 
(normalized) actions for a massles^l scalar field on a square by 

S n = (4> 2 i + a<t>i<t>i+\ - (1 + a)4>i4)i + 2) (2.7) 

i=l,...,4 

one will find that a = — ^ is the only fixed point. This fixed point action agrees with the previous 
fixed point modulo a global rescaling. 

Next we are going to construct an action which has the same invariance holding for more 
complicated boundary data, i.e. where the boundary fields can have a 'kink' at the midpoint 
of the edges. To each midpoint (i, i + 1) between the corners % and i + 1 we introduce the 
'fluctuation field' 7^+1 = 4>ij — \{<pi + <j>i+i)- As initial action we take Hamilton's principal 

5 which we define by an action which evaluated on constant fields gives zero 
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function for the complex contracted out of four squares, denned by taking the fixed point action 
(|2.8p as the action for one square. So we just have to repeat the previous computation keeping 
the eight boundary fields general. 
The resulting action is given by 

5 D 4 = Sq + ^ (0»7i,i+l + <t>ili,i-l - 4>ili+l,i+2 ~ 4>ili-l,i-2) + 
l<i<4 

^ (ai7f )i+ i + 027t,t+l7i+l,i+2 + «27i,i+l7i-l,i + ^«37i,i+l7«+2,i+3) (2-8) 
l<i<4 

with coefficients ap,/3 = 1,2,3, which flow under the coarse graining procedure: we glue four 
of such squares, solve for the fields on the (now five) inner vertices, and determine Hamilton's 
principal function on the subset of boundary data, where eight of the sixteen boundary fields 
are fixed as suitable averages of the eight remaining fields, see figure [2J We take Hamilton's 
principal function evaluated on such data as the new action for the square. 

The fixed point values for the coefficients as, /3 = 1, 2, 3 can be given as the roots of a higher 
order polynomial, the numerical values are given by 

ai : 2.5833 -> 2.2756, a 2 : -1.5 -)■ -1.2003, a 3 : -0.1666 -> -0.3427 (2.9) 

where we gave the values for the initial Hamilton's principal function and the action/ Hamilton 
principal function at the fixed point. 

One will notice that the zeroth order part - the part quadratic in the (j> fields - does not 
change (if one starts with the previously computed fixed point action) , nor does the part coupling 
the (f) and the 7 fields. We want to emphasize that this is a special property of the dynamics 
and choice of coarse boundary data considered here. To appreciate this point assume that we 
would treat the action (|2.8f) in a perturbative expansion with small inhomogenities 7. To this 
end we introduce an expansion parameter e so that the action has the form 

S = \4>-M + N -7+ \<?~i-K- 7 (2.10) 

where M, N, K denote the quadratic forms specifying the dynamics. For the solutions to the 
equation of motions we make the ansatz 

= fy + ety , 7 = °7 (2.11) 

which have then to satisfy the equations of motion 

= M-% + e (M ■ 1 4> + N- ° 7 ) , = eV N + e 2 ■ N + K . (2.12) 

We can see that a perturbative solution is only possible if %■ N = for solutions of the equations 
M ■ <p = 0. This is a consistency conditions involving both M - that is the definition of the 
dynamics at zeroth order - and N - that is the coupling between zeroth order and first order 
variables. Similar consistency conditions occur for the discretization of systems with gauge 
symmetries in situations where these are broken by discretization, see the discussion in [21j . 

If the consistency conditions are satisfied the zeroth order solution will not be affected by 
the higher order terms (i.e. there are no backreaction effects). Remarkably these conditions 
are satisfied for the action (12. 8ft and explain why only the quadratic part in 7 changes. We 
will however see that the first order term will change (slightly) if we include the higher order 
fluctuations. 

We can introduce further fluctuation fields Km+i and «i which are associated to 

vertices subdividing the edges between (i), + 1) and between + 1), (i + 1) respectively. 
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As mentioned previously the first order part, which is quadratic in the 7-fields now does change, 
indicating that the 7 and k variables are not decoupled. The full fixed point action is quite 
lengthy, therefore we will just give the shifted values for the fixed point values of the coefficients 
for the 7 • 7 part 

a x : 2.2756 -> 2.2011 , a 2 : -1.2003 -»• -1.1805 , a 3 : -0.3427 -> -0.3562 (2.13) 

Note that in the action f|2.8j) all fields are coupled to each other, that is the theory is non- 
local if one considers couplings between fields associated to one building block. On the other 
hand we still built the full dynamics out of local actions and building blocks which are glued to 
each other locally, namely by identifying fields on shared boundaries between building blocks. 
There are no couplings between fields on building blocks which do not share boundary fields. 
Such non-local couplings would however arise in a perfect action approach [121 122} HB"]. where 
one keeps the building blocks elementary. This shift of non-locality can be understood as an 
example of the general method to avoid non-local couplings by introducing additional variables. 

It should now be clear how to iterate this procedure. In each iteration we will allow more 
and more general boundary data. In the continuum limit we can hope to obtain Hamilton's 
principal function for the boundary of a square. 

3 Scalar field with potential 

This method is also suitable for non-linear theories, e.g. a massless scalar field with potential 
X(p 4 . Here even the discretization on a triangle does not enjoy any more the invariance under 
subdivisions. Also there are different possibilities to discretize the potential term, for example 
one could take the average of the potentials of the fields or the potential evaluated on the average 
of the fields associated to one building block. We will see that with our method we can find a 
preferred discretization. 

For the action on the square we make the ansatz that the potential term is a general ho- 
mogenous polynomial of fourth degree in the fields. The kinetic term is given by the fixed point 
action for the massless field (|2.6|) : 

Sn,x = S a + 

Xa 2 (&i0f + 6 2 0f &+i + b 2 4?i4>i-i + h<f>i(f>i+i + h$(f>i+2 + b 5 ($(f>i +1 (f>i + 2 + 

i 

h<t>i<j>i-l<t>i-2 + &6<fo-l<^<fo+l + ^b 7 4> 2 (j) 2 +2 + !& 8 0i0 i+ i0i +2 0i+3) (3.1) 

where a is the lattice constant, so that the potential term is weighted with the area of the 
square. We proceed in the same way as in the previous section, that is we determine Hamilton's 
principal function on the subdivided square evaluated on special boundary data, which also here 
we take to be linear along the edges of the square. The only point one has to take care of is 
to replace a — > a/2 for the action of the smaller squares in the composite square in each coarse 
graining step. The full nondinear solution for the field on the inner vertex is a non-polynomial 
function of the boundary fields and in A. To simplify the calculation we can expand in A and 
consider only the first order in A, see section 0] for the second order result. The first order 
calculation implements also a truncation of the flow to polynomial actions up to degree four. 
As the quadratic part is already invariant only the coefficients bp will be changed. The solution 
to the fixed point equations is given by 

h = c, b 2 = c, 63 = c, 64 = \c, 65 = \c, b 6 = fc, 67 = |c, b 8 = |c (3.2) 

where c is a multiplicative constant that just rescales A. (For the initial action (|3.ip with 
61 = 1, bi^\ = we flow to c = 0.16.) 
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We see that we can define in this way a preferred discretization of the potential term for a 
given lattice. In the next section we will discuss whether this discretization encodes already the 
continuum result for edge-wise linear boundary data. 

4 Larger coarse graining steps 

In this section we are going to discuss the dependence of the sequence of fixed point action on 
changing the coarse graining scheme, in particular by taking larger coarse graining steps. This 
is actually very much related to the issue we discussed in section [5J namely whether the lower 
order results changes if we take higher order fluctuations (more complicated boundary data) 
into account. 

We consider the scalar field (with \(f) A potential) on a square with initial action (|2.5p and 
subdivide not only into four but into sixteen squares. Again we put special boundary data. 
For the simplest case, we choose boundary fields which are linear along the edges of the bigger 
square. For the inclusion of the fluctuations 7 we allow a kink - described by the 7-fields - at 
the middle of the edges. 

In general the fixed point actions will change as compared to the case, where we consider 
only the blocking of four squares into one square. The zeroth order result, i.e. the part which 
is quadratic in the fields 0, does not change - neither if we include fluctuations 7 or k into the 
boundary nor if we change the step size of the coarse graining procedure. 

It is instructive to understand why this happens. The reason is in the choice of coarse 
boundary data - as for the case of the triangle in section ([2]) the piecewise linear boundary data 
are induced by a continuum solution which we can naturally obtain by extending the solution 
for the inner vertices and is given by 

cp{x, y) = ^(1 - \x - xi\) (1 - \y - yi\)(t>i . (4.1) 

i 

Here (x, y) are orthogonal coordinates on the larger square normalized such that the coordinate 
distance between neighbouring vertices of the larger square is equal to one. (xi,yi) are coordi- 
nates of the vertex i. This solution holds for the actions (|2.6j) with any value of a and A = 
and for the subdivision into any number of smaller squares. Indeed in the case that we have a 
larger square built out of more than four squares also the solution along the inner edges will be 
a linear function of some lengths parameter along this edge. 

This reproduction of our choice of coarse boundary data along the inner boundaries is re- 
sponsible for the stability of the zeroth order result under inclusions of higher order fluctuations. 
It also means that we actually found a continuum result: the fixed point action coincides 
with the continuum Hamilton's principal function evaluated on edge-wise linear boundary data. 

The exact reproduction of our notion of coarse boundary data on the inner boundaries does 
not hold anymore if we include either the 7 fluctuations or the potential term Xcj) 4 . For the latter 
note that the terms linear in A, which we computed in section ([3|) will again remain stable both 
under the inclusion of 7 fields and under changing the blocking size from four to sixteen squares. 
The reason is that for the computation of Hamilton's function to first order in A one needs only 
the zeroth order (in A) solution. Therefore the change will only start with the second order of A. 
We therefore computed the fixed point action for the scalar field with Xcp 4 potential to second 
order in A, which is equivalent to a truncation to six order polynomials in the fields. We can 
compare the results with and without inclusion of 7 fluctuations into the boundary data, and 
with a blocking of four squares into one or with a blocking of sixteen squares into one. The final 
fixed point actions are very lengthy as for instance at second order all six degree polynomials 
in the eight boundary fields (4>i and 7^) will appear. We will therefore just give the coefficients 
for the ("/ij) n ,n = 2,4,6 terms in the action and the coefficient for the (\>\ term for the various 
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coarse graining schemes. For comparison in all cases the coefficient in front of (ftf is | and the 
coefficient in front of a 2 X(pj is b\ = 0.16. 





4 -> 1; (f>i 


4 ->• 1; <f)i,jij 


16 -> 1; <fn 


16 -> 1; 0j,7y 
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0.001020 
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(4.2) 



As we can see from the \ 2 4>^ coefficients (and the other A 2 coefficients in the <\> fields not 
displayed here) the difference for the fixed point actions which include 7-boundary fields are 
smaller than for the cases without the 7-fields. Note also that the fixed point action computed 
from the coarse graining with 16 — > 1 square but without including 7 fields is different from 
the fixed point action computed in the 4 — > 1 scheme but with the inclusion of the 7 fields. 
(Compare also the 7^ coefficient computed with including the k fields, equation (12.13|) and in 
the 16 — >• 1 scheme.) The difference between the fixed point actions in the 4 — > 1 and in the 
16 — > 1 scheme involving the <p and 7 fields is however smaller compared to the case where 
only the cj> fields are included. One would expect to obtain the continuum Hamilton's principal 
function if we consider larger and larger coarse graining steps. (This is just the requirement 
of having a discretization that gives the continuum result in the continuum limit.) In such a 
set-up we would however need to solve in one step for all the bulk fields at once. 

Increasing the number of boundary fields instead of the blocking step size has the advantage 
that the computation proceeds step-wise. For example for the 16 — > 1 scheme one has to solve 
for nine (cj)) fields on the inner vertices, whereas for the 4 — > 1 scheme and inclusion of 7 fields 
we just have to solve for five fields. Both schemes involve however the same number of boundary 
data (prior to restricting the boundary data to special 'coarse' values). The reason why the 
second scheme involves less inner fields, is that we are using already Hamilton's function for the 
square with more boundary data, in which we already solved for the field on the inner vertex. 
This gives exactly the difference of the four fields, whose solutions we obtained in the previous 
computation. In general we have a quadratic growth (in 2D) in the number of fields to solve 
for in the scheme with larger and larger blocking steps as compared to a linear growth in the 
scheme with including more and more boundary data. 

To get an impression of the convergence rate of the lower order parts of the fixed point action 
under inclusion of more and more boundary data we consider a scalar field with mass in the 
4 — > 1 blocking scheme and compute the fixed point action by taking only fa fields, fa and 7^ 
and finally fa,jij,Kuj fields into account. Although we are just considering a quadratic action 
the full fixed point action will involve arbitrary high (even) powers of the mass. To make the 
computation feasible we truncate to m 4 , that is to second order in the perturbation parameter 
m 2 . (As for the A terms the m 2 4> 2 terms of the action will remain stable under inclusion of higher 
order fields.) For comparison in all cases the coefficient in front of <fi 2 is | and the coefficient in 
front of a 2 m 2 (pf is 0.44444. 
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Figure 3: The tensor network scheme for the coarse graining of a scalar field. The original square is 
rotated by 45°, so that now the boundary fields can be associated to the edges of an unrotated square. 
The gluing of these new squares corresponds to combining tensors of rank four. The corresponding tensor 
network is indicated by blue lines. Bulk fields fa, . . . ,fa are now associated to the inner edges of this 
tensor network. Coarser boundary data correspond to piece wise constant fields. In the figure on the 
right we depicted the set-up for coarse graining with eight boundary fields - fields on pairs of outer legs 
are set to be equal. 



As expected the change in the coefficient of the m 4 a 4 0? term is decreasing, from an order of 
10~ 3 to an order of 10 -4 between the different truncations. Thus we can hope on a convergence 
of the coefficients to some (continuum) value. 

In conclusion we have to expect 'back reaction' contributions from the inclusion of more 
inhomogeneous boundary data onto the terms in the fixed point action involving more homo- 
geneous fields. This applies even for the linear dynamics in the case that the corresponding 
field modes do not completely decouple, as is the case between the 7 and k fields. Such back 
reaction effects could be avoided by changing to decoupled modes, i.e. for a translation invariant 
boundary the Fourier modes. In this case the higher (order) modes will not appear at all in the 
solution for the bulk fields - as long as these modes are set to zero in the boundary data. Hence 
the boundary data reproduce itself on the inner boundaries, as discussed above. 

Of course the construction of decoupled modes will only be possible for free theories. For 
interacting theories we might be able to identify modes, that is a decomposition of boundary 
fields, which only couple weakly with each other. The (coarse) boundary data will only be 
approximately reproduced on inner boundaries, but we can expect a fast convergence of the terms 
involving only lower order modes under inclusion of higher (order) modes into the boundary data. 

The advantage of the scheme proposed here is that one can determine the size of the cor- 
rections at least from one order to the next order, so that one can judge the reliability of the 
approximation scheme. 

5 Other coarse graining schemes 

Here we will discuss a coarse graining scheme with a different choice of coarser boundary data. 
This choice follows more closely the usual understanding of coarse graining, e.g. fields in a region 
are averaged to a new field, constant on the region averaged over. 

Hence we replace for the coarser boundary data 'piecewise linear' fields with 'piecewise 
constant' fields. For this we have also to change the geometrical set-up, as the fields will be now 
associated to the (midpoints of the) edges. Such a set-up corresponds more closely to a tensor 
network renormalization approach, see [23^ l24"t [25] and figure El 

Coarser boundary data are now obtained by setting fields on a subdivided boundary edge to 
be equal. There is one subtlety one has to take care of: to find a fixed point one has to introduce 
a scale s between the new average boundary fields and the to be averaged finer boundary fields, 
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i.e. (pi s = s(4>ii + 4>i2) for the case that the boundary fields agree on pairs of boundary edges 
i2). To find the fixed point action, one has to normalize one of the coefficients in the action 
after each renormalization step to one. At the fixed point one can then determine the scaling 
s by finding the normalization constant. This scale has to be properly taken into account if a 
potential term is added: as we can interpret s as a length scale, fields carry the dimensions of 
inverse lengths. Hence if we have to multiply the terms </> 2 with s -2 , we should multiply terms 
a 4 cp 2 with s 2 whereas terms a 2 (f> 2 are dimensionless and will not be rescaled. 

We determined the fixed point action for the scalar field with mass, in a scheme where four 
squares are blocked into one square. We computed the fixed point action for three different 
set of boundary data: starting with four boundary fields {</>,, % = 1, . . . 4} which we interpret as 
edge-wise constant fields, to a set {(pn, (j>i2, i = 1, • • • ,4} of eight boundary fields which we can 
interpret as piecewise constant fields along the edges, and finally a set {<fo> <fe> ^3, 4>i4\ of sixteen 
boundary fields. We will just display the coefficients for these three cases after restricting the 
boundary values to the simplest case, i.e. setting cpn = fci = ... in the fixed point action. The 
resulting actions have to be normalized again in the way described in the previous paragraph, 
so that the actions can be compared with each other. The initial action is given by 

Su = - + a 2 m 2 <p 2 ) . (5.1) 

i 

In the following table we will display the results for the three different sets of boundary fields 
with the notation 

Sg = y^ j {4>1 + bi(pi4>i + i + \b 2 4)i4>i + 2 + a 2 m 2 (ci4) 2 + c 2 0j0j+i + |c 3 0i0i +2 ) 

i 

+ a i m i {di4> 2 + d 2 0i<fo+i + ^3^1^+2)) • (5.2) 

Because of the rescaling procedure the coefficient of the term <p 2 will be always one, furthermore 
we will have 62 = — 2(1 + b\) which ensures that the massless part vanishes, if evaluated on a 
constant boundary field. 





4 bdry fields 


8 bdry fields 


16 bdry fields 


h 


0.928 


0.945 


0.958 


Cl 


0.493 


0.534 


0.536 


C2 


0.408 


0.370 


0.366 


C3 


0.197 


0.192 


0.191 


di 


-0.189 


-0.203 


-0.261 


d 2 


-0.260 


-0.283 


-0.373 


d 3 


-0.181 


-0.214 


-0.286 



(5.3) 



Prom table ()5.3[) we can see that we do not have stability of the coefficients at zeroth order, 
as was the case for the scheme with piecewise linear boundary fields. We explained in section 
H] that this cannot be expected as piecewise constant fields do not appear as solutions on the 
inner edges. The coefficients for the kinetic part of the action and the ones in front of the m 2 
terms do converge, although slower than in the scheme with piecewise linear boundary fields. 
Moreover the coefficients in front of the m 4 terms show an unstable behaviour. In summary the 
scheme with piecewise linear fields seems to be better suited for the dynamics of a scalar field 
(perturbed around the massless case). 

In general we have to choose carefully a notion of coarse boundary data, which will underly 
the coarse graining scheme. This notion should be adjusted to the dynamics - ideally one should 
capture the relevant degrees of freedom. 
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Another possibility to determine a notion of coarse boundary data is to follow the tensor 
network renormalization algorithm, where the truncation to coarse boundary data is entirely 
determined by the dynamics. This procedure was however developed not for classical theories but 
for quantum ones (with discrete configurations, like Ising spin systems), in which the partition 
function is rewritten as a contraction of a tensor network. Contracting just four tensors T to a 
new tensor T' , this tensor will have the double number of indices, which correspond to the double 
number of boundary data in each coarse graining step in our case. Thus one has to truncate the 
tensor down to its original size as otherwise the number of indices grows exponentially with the 
number of iterations. This truncation is motivated by finding the most relevant modes in the 
tensor T that contribute mostly to the partition function. The general problem in identifying 
these modes is that there is no canonical definition of eigenvalue decomposition for tensors, as 
there is for matrices. Thus a procedure has been designed that introduces auxiliary matrices 
that allow for an eigenvalue decomposition |23t 124] . In general, the identification of the 'most 
relevant' modes leads to field redefinitions, that may also depend on the renormalization step. 
The modes, identified as relevant, are included in the next iteration step, whereas the irrelevant 
ones are discarded. 

As the entire procedure is not so much motivated by finding the boundary wave function 
(see however the discussion in [23]) for given boundary data, one usually does not keep track 
of these field redefinitions. This could however be done in principle. The information gained 
in this way could be used to identify a notion of coarse boundary data that is adjusted to the 
dynamics and therefore might also change during the iteration procedure. The resulting notion 
will of course depend on the truncation scheme for the tensors. On the other hand if one has 
reason to believe that a given notion of coarsening captures the relevant degrees of freedom, this 
notion can be used to define a truncation scheme for the tensor network algorithm. Section [7] 
will discuss a possibility to extend our classical scheme to quantum mechanics, and we will see 
that identifying a proper notion of coarse boundary data becomes even more essential. 

6 Cylindrical consistency 

In the discussion so far a choice of coarse boundary data was essential. We can formalize this no- 
tion using the concept of cylindrical consistency. This is a key idea in loop quantum gravity [2] , 
which allows the construction of a Hilbert space encoding the continuum configurations through 
an inductive limit of a family of Hilbert spaces that are associated to (finite) graphs. There 
is a natural (partial) ordering of graphs into coarser or finer graphs. Cylindrical consistency 
implements the requirement that a computation (of expectation values, transition amplitudes 
etc.) for a given state should not depend on whether this computation is performed in a Hilbert 
space associated to a coarser or finer graph, as long as both graphs are fine enough to describe 
the given state. 

Following this idea we introduce a (partial order) directed set structure on the set B of bound- 
aries b, where b -< b', if b' is a refinement of b. To every boundary we associate a configuration 
space Cb, which typically is a direct product of some basic configuration space C. 

The notion of coarse boundary data will be encoded in the choice of a family of embedding 
maps. For every pair b -< b' we define an injection 



where the image of i^y determines the 'coarse boundary data' in the configuration space Cy ■ 
This family of injections needs to be consistent in the sense that for b -< b' -< b" we should have 



tbb 1 '■ Cb — > Cy 




tbb" = I'Vb" L bb' 



(6.2) 
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Such a consistent family allows the construction of an inductive limit 

Cind = U<V ~ , (6.3) 

which consists of equivalence classes of elements in a disjoint union of configuration spaces C b 
over all boundaries b € £>. Two elements c b and cy are equivalent if there exist a b" with b -< b" 
and b' -< ft" and ibb"{cb) = *w( c fc') holds. In other words two configurations defined on two 
different boundaries are equivalent if there is a finer boundary onto which the two configurations 
can be embedded and if these embeddings happen to agree. 

A family of functions {Fb}beB is cylindrically consistent on the inductive family defined by 
(C b ,i b ), if 

F b = u* w F V: i.e. F b (c) = F v (l w {c)) Vc G C b (6.4) 

where if., is the pullback of i bb >. This just implements that the value of the function F should 
not depend on the representative on which one chooses to evaluate. 

The procedure proposed in this work attempts to construct Hamilton's principal function as 
a cylindrically consistent object: Assume that we do have the continuum Hamilton's function at 
our disposal. In the language we are using here, the continuum is given by the refinement limit 
of the boundaries b, that is on a heuristic level we can imagine the continuum as an element c 
for which b -< c for all finite boundaries in B. Now we can pull-back the continuum Hamilton's 
function to any coarser boundary and the family of function so obtained will be cylindrically 
consistent. To be more precise the continuum limit of a function is actually defined through the 
condition of having a cylindrically consistent family of functions. 

We compute this continuum limit iteratively starting from an action Sb (identified with 
Hamilton's principal function) on some given boundary b. This is used to compute Hamilton's 
principal function on some finer boundary b', we will denote the result by St. For the fixed 
point action we demand that 

L bb'Sb' = &b ■ (6-5) 

If the fixed point action 5£ in (|6.5|) does not depend on the choice of finer boundary ft', then it is 
the ft-component of the cylindrically consistent family of Hamilton's principal functions we are 
looking for: in the limit that b' approaches the continuum we converge to the continuum Hamil- 
ton's principal function and the pull back to the coarser boundary b defines the ft-component of 
the corresponding cylindrically consistent family of functions. 

We have seen that in the case of the massless free scalar field and for the zeroth order 
coefficient such an independence holds, but that in general one can only hope for approximate 
(or a convergence) of results. Here the choice of the finer boundary b' determines the order of the 
approximation for the ft-component of the cylindrical consistent family of Hamilton's functions. 

We should remark that (16. 5[) is to be understood for systems without an explicit scale, such 
as the massless scalar field. Note that gravity, where the metric is a dynamical variable, is in 
this sense also without a scale (if one does not expand on a background) |20j . This also holds 
for parametrized systems, where embedding coordinates are added as dynamical variables [20] . 
In all these cases a notion of scale can be replaced by the notion of coarser and finer boundary 
data. 

In the case of scale dependent actions, which we denote by S bt a, for example the scalar field 
with potential, we have to replace (|6.5[) by 

L bb'^b',La = Sb,La (6-6) 

where the lattice constant in b' is La and Sy a L is Hamilton's principal function computed from 
the action S ba . Even more generally the embeddings could also depend on the two scales 
involved. 
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In the previous scalar field examples, based on piecewise linear fields, the basic configuration 
space is associated to a vertex and given by It. The configuration space for a boundary with 
N vertices is given by M N . We considered the set of boundaries of a square with a given 
number of repeated subdivisions of the edges into two halves, b' is a refinement of b if b' can be 
obtained from b by such a subdivision. Let us consider the case that b' is obtained from b by 
one subdivision. We label the vertices in b by indices i (cyclically ordered) and the additional 
vertices in 6', which are between i and i + 1 (mod the number of vertices in b) by (i, i + 1). The 
injection map is then given by defining the fields on the additional vertices as average of the 
neighbouring vertices, i.e. 



Iterating this procedure one obtains a directed set of boundaries and can define the inductive 
limit of configuration spaces. 

7 Remarks on the quantum theory case 

As mentioned in the previous section the notion of cylindrical consistency is crucial in the 
loop quantization of gravity. We therefore remark that the notion of cylindrical consistency for 
Hamilton's principal function translates into the notion of a cylindrical consistent (path integral) 
measure in the quantum theory, a notion proposed in [9]. Such a cylindrical consistent measure 
would allow to make the path integral well defined [9]. Here we will only sketch the idea and 
point out one of the main issues one would have to address even in the case of free theories. 

Hamilton's principal function associated to a region with a boundary is the classical analogue 
to the path integral over this region. The result of this path integral as a function of the boundary 
data is the (generalized) boundary wave function |28j . Note that - given a prescription for the 
path integral this boundary wave function is the unique physical state associated to this regionlf] 

For the following it is useful if we take a dual point on the path integral and rather see it as 
an (anti-linear) functional on the (kinematical) boundary Hilbert space. That is we assume that 
we can associate to every boundary of a region a kinematical boundary Hilbert space Hbound-, 
that encodes all states that can in principle occur kinematically. A well defined quantum theory 
should then associate an amplitude map to every region, that is it defines a map on every 
boundary Hilbert space: 



As we map to C the co-kernel of this map is one-dimensional, defining the physical wave func- 
tion, which we mentioned before. The maps V for unions of regions have to satisfy consistency 
conditions which are generalizations of the usual composition rule for the path integral [28]. 
Classically these consistency conditions correspond to the additivity of the action (or the com- 
position rule for Hamilton's principal function). 

These considerations also hold for the case that we replace space time with a lattice or 
more general discretization. In this case the boundaries are discrete boundaries b to which we 
associate Hilbert spaces %b- As for the classical case we have to choose injection maps i^y that 
embed the Hilbert space lib into Hyy ■ Such injections can be naturally constructed if we have 
projections iryb '■ Cy — >■ C& on the configuration spaces at our disposal: assuming a polarization 

6 Usually one would think of the path integral to propagate states from one boundary to a second one. Here we 
have only one boundary and the wave function corresponds to the Hartle-Hawking 'no boundary' wave function. 
This is fully consistent with the classical phase space one would associate to such a boundary, as such a classical 
phase space is totally constrained [29].) 



I<bb>(4>l,<p2,<p3,- ■ ■) = (01) +02), 02, \{4>2 +03),03,---) 



(6.7) 
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in which the quantum states are functions on the configuration space we can use the pullback 
of the projections to define injections: 



tbb> '■ Ut — >■ lib' 

i/j b ^ipb' where ^b'ifiv) = ^b^Vbipbf)) ■ (7-2) 

The reader will note that in the classical case we also used injections instead of projections. 
In the case we discussed here we can however obtain projections by 'forgetting' the fluctuation 
variables. I.e. in the case that coarse graining was based on piecewise linear fields and for 
a refinement which doubled the number of fields, {4>i} to {(^,0^+1}, we perform a variable 
transformation 7^ = faj — \ {4>i + 4>j) and define the projection as 

7T&'&(01,712,02,---) = (01,02, •••) • (7.3) 

Such a projection satisfies n^b o tfy ss = Id& where here the t^ lass is the classical embedding map 
for the configuration spaces. 

Such inductions which are based on projections are also used for the construction of the 
Hilbert space in Loop Quantum Gravity [7J, see also |30j . 

A measure (which here will be just understood as linear functional) on a family of inductive 
Hilbert spaces can be defined by providing a representation on each of the Hilbert spaces Tib- 
Such a family of measures {[i}b is cylindrically consistent if 

fJ-bi^b) = Hv^bviipb)) ■ (7.4) 

The Ashtekar-Lewandowki measure [7J [2] for the kinematical Hilbert spaces in loop quantum 
gravity satisfies such consistency relations. 

One way to define the dynamics (and the physical inner product) would be via the path 
integrals acting as functionals on these kinematical Hilbert spaced^ as in (17. ip . Cylindrical 
consistent dynamics would then mean to have a cylindrical consistent family of (anti-linear) 
maps 

V b :M b ^€ . (7.5) 

As an example we can again consider the massless free scalar field. The kinematical Hilbert 
space associated to the boundary with N vertices will be L 2 (R N ), here we will have the case 
N = 4 and N = 8. We define the path integral map for the boundary b of the basic square as^| 

Wb) = W f II d0;e^)^(0i) (7.6) 
NbJ i=w 

where Nf, is a normalization factor. The path integral map for a refinement b' that subdivides 
a square into four would be given by 

Vb>(M = 4; / d 0o II d 0* d 7i,i+i e^'^'^+^^^v+i) (7-7) 

b ' ^ i=l,. ..,4 

where 0o is the scalar field on the inner vertex, "fa+x are the variables introduced in (|7.3p and 
Sy is here the action for four squares obtained by summing the basic action Sb over the four 



technically on some dense subspace of these Hilbert spaces 

8 The complex conjugation of the kinematical wave function can be interpreted as the kinematical wave function 
being associated with the other orientation of the boundary [28]. 
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squares. If we define the injections i^y as pullback of the projections (|7.3p we obtain for the 
pullback of the path integral map from b' to b 

{iwfVv {^)=Vv{Lw{iph)) = ^ [d(j> [] d0id 7M+ i e iS >'^"+^M&) ■ (7.8) 

b ' i=l,...,4 

For cylindrical consistency of the path integral maps defined in this way we would need 

d& e iS ^M4>i) = fdcfr J] d ^ d 7M+l e^^M-fiA)^.) (79) 
6 ' ^ i=l,...,4 

for a suitable class of test states ifo. As we are considering free field theory the integral over the 
7 fields and (fro is a (analytically continued) Gaussian integration. Indeed, as pre-factors arising 
from these integrations can be absorbed into the normalizations iV&, Ny we can just consider the 
extremization of Sy with respect to (fro and then with respect to the fields ji i+i- 

Compared to the classical coarse graining procedure we therefore do not set the 7^+1 = 
but have to find the extrema of Hamilton's principal function H(<pi,jij) = S&'^j, 7j,-, 0q o1 ) with 
respect to these variables 7. We then obtain a new action, by evaluating Hamilton's function 
on these extrema. This new action will depend just on the four fields (fri and we can repeat the 
procedure and in this way obtain a renormalization flow of actions. 

The resulting flow does however not converge to a fixed point. Rather a rescaling fixed 
point is reached, in which the action is rescaled by a given number (smaller than one). Indeed 
compared to the classical case we do not fix the boundary values of 7, but do minimize further 
the action with respect to these variables. The values of the boundary fields 7 obtained in this 
way are however not reproduced, if we put the resulting squares again together, so that we have 
now a variation of the fields 7 in the bulk. This is different from the classical case, where the 
boundary values 7 = are reproduced in a larger square. Indeed if we include in this procedure 
higher order fields the (rescaling) fixed point, which differs from the classical one, will change 
towards the location of the classical fixed point. 

Although one could deal with such a rescaling fixed point, the issue seems rather that we do 
not reproduce the classical result in the semi-classical limit. The key point responsible is the 
choice of embedding map, which also describes the choice of boundary conditions (for the vertices 
carrying 7 fields) imposed on the quantum mechanical wave function, which we are trying to 
find. This point constituted also the main motivation for the introduction of the density matrix 
renormalization group procedure [3T]: find the ground state of a given Hamiltonian in a given 
block with some fixed boundary conditions (for instance that the wave function vanishes on 
the boundary) . Putting the blocks together and finding the ground state again for the larger 
block it might happen that the wave function evaluated on the now inner boundaries differs 
very much from the values prescribed by the conditions on the smaller blocks. Even for the free 
particle in one dimension one will find a ground state wave function that vanishes on the ends 
of the smaller intervals, versus a maximum of the ground state wave function assumed on the 
midpoint of the doubled interval. This issue is solved by a kind of bootstrap approach - also 
the appropriate boundary conditions have to be constructed by considering not only the block 
itself but the block inside an environment (built from the same block). 

For our case (of a free field) we can attempt to reproduce the classical coarse graining by 
choosing an appropriate embedding map. In this case we have to pre-guess the correct vacuum 
for the finer degrees of freedom. We can choose the embedding as 

ifrb{(fri) ^ 4> b {cfri)ex.p(iT('y ■ ^,7 • 7)) . (7.10) 
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The condition (17.91) turns into 



i=l,...,4 " " i=l,. ..,4 

(7.11) 



We can fix T such that the term T{j-<j), 7-7) exactly cancels the 7 dependence in the action Sy. If 
Sb is the (classical) fixed point action, we reproduce also cylindrical consistency in the quantum 
case. (This will lead to a divergence in the integral over 7. In the case of loop quantum gravity 
the corresponding integration is over a compact space. Also here we could (Bohr) compactify R, 
equivalent with a regularization of the integral to a finite interval [— C, C] and a normalization 
by 1/2C, letting then C -> 00.) 

We can interpret the choice of embedding (|7.10p as selecting the correct vacuum for the finer 
degrees of freedom 7. Only test states Vfc 1 e * T which have already the correct semi-classical 
behaviour with respect to the 7 degrees of freedom, are propagated in the same way on the finer 
boundary, as the states tpf, on the coarser boundary. 

In conclusion the choice of embedding seems to be even more essential in the quantum 
case. For the construction of these embedding maps methods from the density matrix renor- 
malization group or tensor network based renormalization [23} I2"4"[ [32] , which can be seen as a 
descendent of the density matrix renormalization group procedure, could be helpful. We would 
not have encountered problems with (17. 9p , if the finer and coarser degrees of freedom where 
completely decoupled. Achieving such a decoupling is one of the motivations behind the ap- 
proach of entanglement renormalization [32] . The tensor networks can in this case be understood 
as parametrizing an ansatz for a physical wave function. 

Here we started with an embedding map inspired from loop quantum gravity, which we ar- 
gued is not an appropriate one for coarse graining the free scalar field (on a flat background). 
This does not mean that the loop quantum gravity embedding map might not be useful for 
coarse graining - as we argued the embedding should depend on the system and the dynamics 
it describes. Indeed there are uniqueness results available for the loop quantum gravity embed- 
ding, and the vacuum underlying this embedding [33], which are based on the requirement of 
a (spatially diffeomorphism symmetry covariant) irreducible representation of the kinematical 
observable algebra and a (spatially) diffeomorphism invariant vacuum. On the other hand, it is 
not guaranteed that this embedding will allow for the construction of a cylindrically consistent 
measure encoding the dynamics via a path integral, in the way we outlined here, i.e. via a coarse 
graining procedure. As we have seen in the scalar field example we might need already some 
information on the physical vacuum for the finer degrees of freedom. The question is whether 
we can encode this physical vacuum into an inductive limit structure for the physical Hilbert 
space. The crucial question will be whether such a structure exists and how it relates to the 
inductive limit structure of the kinematical Hilbert space. Note also that alternatives for the 
kinematical Hilbert space have been explored |34[ [30] . 



8 Summary, open issues and outlook 

We proposed a coarse graining scheme, in which the central object is Hamilton's principal 
function for the classical case and the physical wave function for the quantum case. Both can be 
taken to define an improved or perfect discretization |12 [ I14 [ [T6]: Hamilton's principal function 
can be used as the action associated to a building block with the corresponding boundary and 
the physical wave function can be used as the path integral amplitude for this building block. 
The main properties for this coarse graining scheme are: 

9 In gravity physical states have to be annihilated by the diffeomorphism and Hamiltonian constraints 
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• Following the tensor network renormalization [23], instead of introducing non-local cou- 
plings between building blocks we introduce building blocks with more and more boundary 
data. The non-localities are shifted to occur in-between the boundary data associated to 
one building block. An advantage of this scheme is that using these building blocks to 
define a discretization the amount of non-locality (i.e. the maximal lattice distance at 
which fields are coupled) can be chosen. This choice also determines the truncation of the 
renormalization flow, as couplings over larger lattice distances are neglected. 

• The quality of this truncation can be tested by considering the higher order corrections. 
The approximations to Hamilton's function for a given boundary b are specified by finer 
boundaries b', see formula (|6,5p . Such an approximation can be understood as an effective 
action which takes the corrections into account, that arise from the additional degrees of 
freedom in b' as compared to b. Taking the refinement limit of b' one obtains the continuum 
limit of Hamilton's function for boundary data described by b. 

• Crucial for the coarse graining scheme is the choice of coarse boundary data encoded into 
a family of embedding maps. These embedding maps provide also an interpretation of 
discrete boundary data as configurations in a continuum theory. The embedding maps 
have to be adjusted to the dynamics, ideally the notion of coarse boundary data should 
be reproduced by the dynamics of the system. Equivalently the choice of coarse boundary 
data should correspond to an ordering of the variables into finer and coarser degrees of 
freedom which are dynamically only weakly coupled. 

• The scheme proposed here seems to be especially suited for systems without an explicit 
(lattice) scale, which includes gravity systems [20]. Here a notion of scale can be replaced 
by a notion of coarser and finer boundary data. A certain truncation might give reliable 
result for boundary data up to a given complexity. 

There are many open issues left to be understood, in particular for the quantum theory 
case. The main point is how to select best the embedding maps. Here methods developed in 
condensed matter and quantum information theory [31[ [32] will be helpful. One question to 
address will be in which situations it is appropriate to determine the truncation scheme (or 
embedding maps) completely by the dynamics, as is done in [23[ [2^1 [25] [32] or whether some 
predefined embedding might be sufficient. For the investigation of phase transition one would 
rather expect to need truncations that are determined by the dynamics. The second possibility 
might arise if we aim at determining a perfect discretization of a system in a well defined phase. 
This will be investigated for simplified models |35[ 125] in [36] , 

For instance the Ashtekar-Lewandowski measure (and the embedding maps and kinematical 
vacuum on which this measure is based), used in Loop Quantum Gravity, corresponds to the 
high temperature fixed point of lattice gauge theory [9] . Here an interesting question is whether 
we could construct a Hilbert space and measure corresponding to the low temperature fixed 
point of lattice gauge theory, which is actually given by a topological theory, namely BF- 
theory, see for instance [35] ■ BF theory underlies spin foam quantization, hence it would be 
natural to consider a coarse graining procedure with corresponding embedding maps. In this 
case the ('kinematical') vacuum state should correspond to a physical state of BF theory, which 
in (2+1) D even coincides with the physical state of gravity. As shown in [37] the construction 
of this physical state can be obtained starting from the usual loop quantum gravity Hilbert 
space. Hence we see that it is actually possible to obtain a physical vacuum, starting from a 
kinematical Hilbert space, which is based on embeddings which rather encode a kinematical 
vacuum describing a very different phase. The question is however, which are the appropriate 
embedding maps to choose for the coarse graining procedure proposed here. For other LQG 
related Hilbert spaces, based on 'condensate' vacua, describing non-degenerate geometries, see 
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Another important question to address is whether one can estimate and obtain a bound for 
the 'back reaction' effects, i.e. the change of the lower order results, obtained for very coarse 
boundary data, if finer and finer boundary data are included. This would already entail the 
existence of such a continuum limit, which in itself is a highly non-trivial question, in particular 
if it comes to the quantum theory. The advantage of the scheme here is that it provides a 
systematic way to improve the truncation in steps, and that the higher order corrections can be 
checked explicitly. 

Let us mention some possible applications of this scheme. 

The main motivation is the construction of perfect discretizations for (quantum) gravity 
models. Perfect discretizations have been proposed for gravity [131 E] in order to avoid the 
breaking of diffeomorphism symmetry for generic discretizations |19j . The notion of diffeomor- 
phism symmetry for discretizations [27\ [T9] is a very powerful one as it implies discretization 
(or triangulation) independence \17 \ \20 \ I21|. In the scheme proposed here the requirement 
of discretization or triangulation independence is imposed by computing fixed point actions at 
every step. We can then hope to obtain triangulation independence up to a given measure of 
fineness of the boundary data. In its simplest form this is already realized for (classical) Regge 
gravity [6], namely for boundary data that lead to flat |38U17| or constantly curved solutions [39] 
(for gravity with a cosmological constant). Correspondingly a natural choice of an embedding 
map would be to implement piecewise flat or piecewise constantly curved geometries as coarser 
boundary configurations, see also [49]. 

More generally the methods proposed here could be used to define improved discretizations 
for partial differential equations, if these can be derived from an action principle. The use of 
renormalization methods to obtain improved discretizations has been proposed for instance in 
|40[I41|. in particular in |41| there are also embedding (and projection) maps used. The difference 
with the scheme here is the utilization of the action principle, the emphasis on interpreting 
discrete boundary data as elements in finer configuration spaces, as well as the possibility to 
improve the discretization in steps and to control the amount of non-local couplings. 

A particular problem for discretizations, related to the breaking of diffeomorphism symmetry 
is the loss of energy momentum conservation, see for instance \42\ I21j. Perfect discretizations 
would avoid such a breaking as these are defined to display continuum physics. This problem 
can be studied in parametrized field theory, where diffeomorphism symmetry is added through 
the introduction of embedding variables |43j . The problem of energy conservation translates 
then into the problem of broken diffeomorphisms and in the canonical set-up to an anomalous 
constraint algebra. With the methods proposed here we can hope to obtain a discretization that 
is perfect for a certain class of boundary data and for which anomalies in the constraint algebra 
can be avoided. On a general method to obtain a discrete canonical formalism with the same 
(diffeomorphism) symmetries as the covariant one, see |19 [ I44 [ [29]. in particular [29] is applicable 
to boundaries with varying number of variables, which are also occurring here. 

Of course we hope that the scheme proposed here can also be applied to quantum gravity 
models, in particular spin foams [U [5]. Although spin foams, being seen as a quantum version 
of Regge calculus have been originally defined on simplicial lattices only, recently models on 
arbitrary lattices (two-complexes) have been constructed [45, 46J. These also provide an example 
of 'more complicated' (non-simplicial) building blocks, which allow the allocation of finer and 
finer boundary data. We propose here that the amplitudes for such building blocks can be 
constructed through a coarse graining process as outlined in this work. See also [25] for a 
formulation of spin foams as tensor networks and a coarse graining procedure for simplified 
models. There are also formulations available which define spin foams as generalized lattice gauge 
theories |47[ l4~8l [91 06], which again can be extended to any two-complex. These formulations 
might allow to systematically investigate which notion of coarse boundary data, or embedding 
maps, will lead to interesting results. In particular [35[ [25l |4"6] provide a range of simplified 
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models for which this question can be more easily explored. 

The standard LQG embedding provides a well defined choice and is strongly motivated by 
(spatial) diffeomorphism symm etry . As outlined in [39], at least on the classical level there are 
different interpretations possible 10 ! of how the discrete data are to be understood in a continuum 
phase space. We suggest here that an embedding of discrete data into a continuum configuration 
or Hilbert space should be determined by the dynamics. Such an embedding would be rather 
associated to a physical vacuum and a physical Hilbert space - as we indicated in section [7] it 
would encode information on the vacuum wave function for the finer degrees of freedom. 

Finally, we would like to mention that the notion of coarse and finer boundary degrees 
of freedom can be brought together with symmetry reduction and applications in (quantum) 
cosmology. Very coarse states would correspond to homogeneous configurations, finer states 
would describe inhomogenities. Indeed in the context of loop quantum cosmology [50] there 
are different suggestions how the kinematics and dynamics can be obtained from the full theory 
[34 1 I51 4 I52] and how to include corrections from inhomogenities, with |52j discussing (cylindrical) 
consistency conditions for the truncation of a given theory to homogeneous degrees of freedom. 

Acknowledgements 

I am very much indebted to Benjamin Bahr for many discussions, comments on the draft and 
suggesting the notion of cylindrical consistency for the coarse graining scheme. I furthermore 
thank Jan Ambj0rn, Laurent Freidel, Frank Hellmann, Carlo Rovelli, Lee Smolin and Sebastian 
Steinhaus for discussions. Research at Perimeter Institute is supported by the Government 
of Canada through Industry Canada and by the Province of Ontario through the Ministry of 
Research and Innovation. 



References 

[1] R. Loll, "Discrete Approaches to Quantum Gravity in Four Dimensions", Living Reviews in Relativity 13 

(1998) online article: http://relativity.livingreviews.org/Articles/lrr-1998-13/ . 
[2] T. Thiemann, "Modern Canonical Quantum General Relativity", (Cambridge University Press, Cambridge 

2005), |gr-qc/01 10034] , 

[3] B. Bahr, R. Gambini and J. Pullin, "Discretisations, constraints and diffeomorphisms in quantum gravity," 
SIGMA 8 (2012) 002 [arXiv: 1111. 18791 [gr-a c] ] . 

M. Dupuis, J. P. Ryan and S. Speziale, "Discrete gravity models and Loop Quantum Gravity: a short 

review, " arXiv: 1204.5394 [gr-qc]. 
[4] C. Rovelli, "Quantum Gravity" (Cambridge University Press, Cambridge 2004). 
[5] A. Perez, "The Spin Foam Approach to Quantum Gravity," arXiv:1205.2019 [gr-qc]. 
[6] T. Regge, "General relativity without coordinates," Nuovo Cim. 19 (1961) 558. 

T. Regge and R. M . Williams, "Discrete structures in gravity," J. Math. Phys. 41 (2000) 3964 

|arXiv:gr-qc/0012035] . 

[7] A. Ashtekar and C. J. Isham, "Representations of the holonomy algebras of gravity and nonAbelian gauge 
theories," Class. Quant. Grav. 9 (1992) 1433 [hep-th/9202053 



A. Ashtekar and J. Lewandowski, "Representation theory of analytic holonomy C* algebras," gr-qc/9311010 
A. Ashtekar and J. Lewandowski, "Projective techniques and functional integration for gauge theories," J 
Math. Phys. 36 (1995) 2170 [gr-qc/9411046] . 

A. Ashtekar and J. Lewandowski, "Differential geometry on the space of connections via graphs and projective 
limits," J. Geom. Phys. 17 (1995) 191 lhep-th/9412073] . 

A. Baratin, B. Dittrich, D. Oriti and J. Tambornino, "JN on-commutative flux representation for loop quantum 
gravity," Class. Quant. Grav. 28 (2011) 175011 [arXiv: 1004.34501 [hep-th]]. 



10r 



These different interpretations arise trough an interpretation of the discrete phase space as a reduction 
of a continuum phase space by a certain set of gauge symmetries. Different gauge fixings lead to different 
interpretations. 



20 



T. Thiemann, "Anomaly- free formulation of non-perturbative, four-dimensional Lorentzian quantum grav- 
ity," Phys. Lett. B 3 80 (1996) 257 [arXiv:gr-qc/9606088] . "Quantum spin dynamics (QSD)," Class. Quant. 
Grav. 15 (1998) 839 |arXiv:gr-qc/9606089j . 

B. Bahr, "Operator Spin Foams: holonomy formulation and coarse graining," J. Phys. Conf. Ser. 360 (2012) 
012042 |arXiv:1112.3567 [gr-qc]]. 

J. A. Zapata, "Loop quantization from a lattice gauge theory perspective," Class. Quant. Grav. 21 (2004) 
L115 [gr-qc/0401109] . 

E. Manrique, R. Oeckl, A. Weber and J. A. Zapata, "Loop quantization as a continuum limit," Class. Quant. 
Grav. 23 (2006) 3393 [hep-th/05 11222] . 

J. Martinez, C. Meneses and J. A. Zapata, "Geometry of C-fiat connections, coarse graining and the contin- 
uum limit," J. Math. Phys. 46 (2005) 102301 |hep-th/0507039] . 

C. Rovelli, "On the structure of a background independent quantum theory: Hamilton function, transition 
amplitudes, classical limit and continuous limit," arXiv:1108.0832 [gr-qc]. 

P. Hasenfratz and F. Niedermayer, "Perfect Lattice Action For Asymptotically Free Theories," Nucl. Phys. 
B 414 (1994) 785 [arXiv:hep-lat 79308004] . 

P. Hasenfratz, "Prospects for perfect actions," Nucl. Phys. Proc. Suppl. 63 (1998) 53 [arXiv:hep-lat/9709110] . 
K. Giesel and T. Thiemann, "Algebraic Quantum Gravity (AQG). I. Conceptual Setup," Class. Quant. Grav. 
24 (2007) 2465 [gr-qc/0607099] . 

B. Bahr and B. Dittrich, "Improved and Perfect Actions in Discrete Gravity," Phys. Rev. D 80 (2009) 124030 
|arXiv:0907.4323l [gr-qc]]. 

B. Bahr, B. Dittrich and S. He, "Coarse graining free theories with gauge symmetries: the linearized case," 
New J. Phys. 13 (2011) 045009 [arXiv: 101 1.36671 [gr-qc]]. 

B. Bahr, B. Dittrich and S. Steinhaus, "Perfect discretization of reparametrization invariant path integrals," 
Phys. Rev. D 83 (2011) 105026 [arXiv: 110 1.47751 [gr-qc]]. 

B. Dittrich and S. Steinhaus, "Path integral measure and triangulation independence in discrete gravity," 
Phys. Rev. D 85 (2012) 044032 [arXiv: 11 10.68661 [gr-qc]]. 

B. Dittrich, "Diffeomorphism symmetry in quantum gravity models," Adv. Sci. Lett. 2 (2009) 151, 
|arXiv:0810.3"594l [gr-qc] .] 

B. Bahr and B. Dittrich, "(Broken) Gauge Symmetries and Constraints in Regge Calculus," Class. Quant. 
Grav. 26 (2009) 225011 [arXiv:0905. 16701 [gr-qc]]. 

B. Bahr and B. Dittrich, "Breaking and restoring of diffeomorphism symmetry in discrete gravity," THE 
PLANCK SCALE: Proceedings of the XXV Max Born Symposium, Wroclaw (Poland), 29 June - 3 July, 
2009, edited by J. Kowalski-Glikman et. al. ,pp. 10-17, [arXiv:0909.5688l [gr-qc].] 

C. Rovelli, "Discretizing parametrized systems: the magic of Ditt-invariance," arXiv:1107.2310 [hep-lat]. 
B. Dittrich, "How to construct diffeomorphism symmetry on the lattice," arXiv:1201.3840 [gr-qc]. 

W. Bietenholz, "Perfect scalars on the lattice," Int. J. Mod. Phys. A 15 (2000) 3341 [arXiv:hep-lat/9911015] . 
M. Levin and C. P. Nave, "Tensor Renormalization Group Approach to Two-Dimensional Classical Lattice 
Models" Phys. Rev. Lett. 99 (2007) 120601 [arXiv:cond-mat/0611687] . 

Z. -C. Gu and X. -G. Wen, "Tensor-Entanglement-Filtering Renormalization Approach and Symmetry Pro- 
tected Topological Order," Phys. Rev. B 80 (2009) 155131 [arXiv:0903.T069l [cond-mat.str-el]]. 
B. Dittrich, F. C. Eckert and M. Martin-Benito, "Coarse graining methods for spin net and spin foam 
models," New J. Phys. 14 (2012) 035008 [arXiv:1109.4927l [gr-qc]]. 

B. Dittrich and F. C. Eckert, "Towards computational insights into the large-scale structure of spin foams," 
J. Phys. Conf. Ser. 360 (2012) 012004 [arXiv:llll. 09671 [gr-qc]]. 

M. Wardetzky, S. Mathur, F. Kalberer, E. Grinspun, "Discrete Laplace operators: No free lunch," Proc. 
Symp. Geom. Proc. 2007 

M. Rocek and R. M. Williams, "Quantum Regge Calculus," Phys. Lett. B 104 (1981) 31. 

R. Oeckl, "A 'General boundary' formulation for quantum mechanics and quantum gravity," Phys. Lett. B 



575 (2003) 318 [hep-th / 306025 . F. Conrady and C. Rovelli, "Generalized Schrodinger equation in Euclidean 
field theory," Int. J. Mod. Phys. A 19 (2004) 4037 [hep-th/0310246] . 

B. Dittrich and P. A. Holm, "Canonical simplicial gravity," Class. Quant. Grav. 29 (2012) 115009 
|arXiv:1108.i974l [gr-qc]]. 

B. Dittrich, C. Guedes, D. Oriti, "On the space of generalized fluxes for loop quantum gravity", to appear 
S. R. White, "Density matrix formulation for quantum renormalization groups" Phys. Rev. Lett. 96 (1992) 
2863. 

U. Schollwock, "The density-matrix renormalization group" Rev. Mod. Phys. 77 92005) 259 
[32] G. Vidal, "Entanglement renormalization," Phys. Rev. Lett. 99, (2007), 220405 

G. Evenbly, G. Vidal, "A real space decoupling transformation for quantum many body systems," arXiv: 
1205.0639 



21 



[33] C. Fleischhack, "Representations of the Weyl Algebra in Quantum Geometry," to appear in Commun. 
Math. Phys. IarXiv:math-ph/0407006 J. Lewandowski, A. Okolow, H. Sahlmann and T. Thiemann, "Unique- 
ness of diffeomorphism invariant states on holonomy-flux algebras," Commun. Math. Phys. 267 (2006) 703 
|arXiv:gr-qc/0504147] . 

[34] T. A. Koslowski, "Dynamical Quantum Geometry (DQG Programme)," arXiv:0709.3465 [gr-qc]. 

T. Koslowski and H. Sahlmann, "Loop quantum gravity vacuum with nondegenerate geometry," SIGMA 8 

(2012) 026 |arXiv: 1109.46881 [gr-qc]]. 
[35] B. Bahr, B. Dittrich and J. P. Ryan, "Spin foam models with finite groups," arXiv:1103.6264 [gr-qc]. 
[36] B. Dittrich, S. Steinhaus, "Coarse graining spin net models with different vacua," to appear 
[37] K. Noui and A. Perez, "Three dimensional loop quantum gravity: Physical scalar product and spin foam 

models," Class. Quant. Grav. 22 (2005) 1739 [arXiv:gr-qc/0402110] . 
[38] B. Dittrich and J. P. Ryan, "Phase space descriptions for simplicial 4d geometries," Class. Quant. Grav. 28 

(2011) 065006 |arXiv:0807.2"806 1 [gr-qc]] . 
[39] B. Bahr and B. Dittrich, "Regge calculus from a new angle," New J. Phys. 12 (2010) 033010 [arXiv:0907.4325l 

[gr-qc]]. 

[40] N. D. Goldenfeld, A. McKane, Q. Hou, "Block spins for partial differential equations", J. Stat. Phys. 93 
(1998), 699 

[41] A. Degenhard, J. Rodriguez-Laguna, "Towards the evaluation of the relevant degrees of freedom in non-linear 
partial differential equations, J. Stat. Phys. 106 (2001), 1093. 

A. Degenhard, J. Rodriguez-Laguna, "Real-space renormalization-group approach to field evolution equa- 
tions, Phys. Rev. E65 (2001) 036703 

[42] J. Marsden and M. West "Discrete mechanics and variational integrators" Acta Numerica 10, 357 (2001) 
[43] C. J. Isham and K. V. Kuchar, "Representations Of Space-time Diffeomorphisms. 1. Canonical Parametrized 

Field Theories," Annals Phys. 164 (1985) 288. 
[44] B. Dittrich and P. A. Hohn, "From covariant to canonical formulations of discrete gravity," Class. Quant. 

Grav. 27 (2010) 155001 |arXiv:0912.18T7 [gr-qc]]. 
[45] W. Kaminski, M. Kisielowski and J. Lewandowski, "Spin-Foams for All Loop Quantum Gravity," Class. 

Quant. Grav. 27 (2010) 095006 [Erratum-ibid. 29 (2012) 049502] |arXiv:0909.0939l [gr-qc]]. 

B. Bahr, F. Hellmann, W. Kaminski, M. Kisielowski and J. Lewandowski, Class. Quant. Grav. 28 (2011) 
105003 |arXiv:1010.4"787l [gr-qc]]. 

[46] B. Bahr, B. Dittrich, F. Hellmann, W. Kaminski, "Holonomy spin foam models: definition and coarse 
graining", to appear 

B. Bahr, B. Dittrich, F. Hellmann, W. Kaminski, "Holonomy spin foam models: Universal boundary Hilbert 
spaces", to appear 

[47] H. Pfeiffer, "Dual variables and a connection picture for the Euclidean Barrett-Crane model," Class. Quant. 



Grav. 19 (2002) 1109 gr-qc/0112002 



[48] E. Magliaro and C. Perini, "Local spin foams," arXiv: 1010.5227 [gr-qc]. 

[49] L. Freidel, M. Geiller and J. Ziprick, "Continuous formulation of the Loop Quantum Gravity phase space," 

larXiv:1110.4"833l [gr-qc] . 
[50] M. Bojowald, "Loop quantum cosmology," Living Rev. Rel. 11 (2008) 4. 

M. Bojowald, "Loop Quantum Gravity and Cosmology: A dynamical introduction," arXiv:1101.5592 [gr-qc]. 

A. Ashtekar and P. Singh, "Loop Quantum Cosmology: A Status Report," Class. Quant. Grav. 28 (2011) 
213001 [arXiv: 1108.08931 [gr-oc]]. 

[51] M. Bojowald and H. A. Kastrup, "Quantum symmetry reduction for diffeomorphism invariant theories of 
connections," Class. Quant. Grav. 17 (2000) 3009 |hep-th/9907042| . 

J. Brunnemann and T. Thiemann, "On (cosmological) singularity avoidance in loop quantum gravity," Class. 

Quant. Grav. 23 (2006) 1395 [gr-qc/05050 32]. 

J. Brunnemann and C. Fleischhack, arXiv:0709.1621 [math-ph]. 

B. Dittrich and J. Tambornino, "Gauge invariant perturbations around symmetry reduced sectors of gen- 
eral relativity: Applications to cosmology," Class. Quant. Grav. 24 (2007) 4543 [gr-qc/0702093| [GR-QC]]. 
J. Brunnemann and T. A. Koslowski, "Symmetry Reduction of Loop Quantum Gravity," Class. Quant. Grav. 
28 (2011) 245014 |arXiv:1012.0"053l [gr-qc]]. 

C. Rovelli and F. Vidotto, "Stepping out of Homogeneity in Loop Quantum Cosmology," Class. Quant. Grav. 
25 (2008) 225024 |arXiv:0805.4585l [gr-qc]] . 

E. Bianchi, C. Rovelli and F. Vidotto, "Towards Spinfoam Cosmology," Phys. Rev. D 82 (2010) 084035 
|arXiv:1003.3"483l [gr-qc]] . 

E. F. Borja, J. Diaz-Polo, I. Garay and E. R. Livine, "Dynamics for a 2-vertex Quantum Gravity Model," 
Class. Quant. Grav. 27 (2010) 235010 |arXiv: 1006.24511 [gr-qc]] . 

E. R. Livine and M. Martin-Benito, "Classical Setting and Effective Dynamics for Spinfoam Cosmology," 
larXiv:1111.2867l [gr-qc] . 

[52] F. Hellmann, "On the Expansions in Spin Foam Cosmology," Phys. Rev. D 84 (2011) 103516 |arXiv:1105.T334l 
[gr-qc]]. 



22 



